FTUV/99-52, IFIC/99-54 

N— quantum approach to quantum field theory at finite T and //: 

the NJL model 



Yasuo Umino and Vicente Vento 
Instituto de Fisica Corpuscular - C.S.I.C. 
Departamento de Fisica Tedrica, Universitat de Valencia 
E-46IOO Burjassot, Valencia, Spain 
(February 1, 2008) 



Abstract 

We extend the N-quantum approach to quantum field theory to finite temper- 
ature (T) and chemical potential (//) and apply it to the NJL model. In this 
approach the Heisenberg fields are expressed using the Haag expansion while 
temperature and chemical potential are introduced simultaneously through 
a generalized thermal Bogoliubov transformation. Known mean field results 
are recovered using only the first term in the Haag expansion. In addition, 
we find that at finite T and in the broken symmetry phase of the model the 
mean field approximation can not diagonalize the Hamiltonian. Inclusion of 
scalar and axial vector diquark channels in the SU(2)f g) SU(3) C version of the 
model can lead to a lowering of the vacuum energy density. We discuss how 
to go beyond the mean field approximation by including higher order terms 
in the Haag expansion. 
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I. INTRODUCTION 



The N-quantum approach (NQA) to quantum field theory developed by O.W. Greenberg 
0II is a method to solve the operator equations of motion by expanding the Heisenberg fields 
as a complete set of normal-ordered and on-shell asymptotic (in- or out-) fields, including 
the asymptotic fields for bound states. This expansion is known as the Haag expansion || . 
Although the NQA is yet to be applied to gauge theories successful applications thus far 
indicate that the method can be very powerful in some cases (see for a summary of the 
applications and references). In particular NQA is well suited to describe bound states. 

Greenberg and collaborators applied the NQA to study dynamical symmetry breaking 
involving the condensates of composite fields in the NJL model [|J both without || and with 
U isospin. Their method differed from the conventional approach to the NJL model in 
several distinct ways: 

• The NQA uses retarded amplitudes (for a Haag expansion using in-fields) with all but 
one leg on-shell instead of off-shell Feynman amplitudes. 

• It allows one to work directly with the fields and states of broken symmetry without 
any reference to those of the unbroken-symmetry vacuum. 

• The explicit structure of the vacuum need not be specified. 

• Bound state amplitudes are calculated directly rather than searching for bound state 
poles in scattering amplitudes. 

In both H and || it was shown that the well known gap equation for the NJL model 
can be obtained by using only the first term in the Haag expansion in the operator equation 
of motion. Analytical expressions for amplitudes and masses for various bound states were 
found by going one order higher in the expansion and the resulting masses agreed with the 
original results |4| except for the mass of the vector bound state. The form factor of the 
Nambu-Goldstone boson was also examined by extending the Haag expansion to third order 
in the one-loop approximation ||. 

Since the NJL model has been used extensively to study the consequences of chiral 
symmetry breaking in hot and dense matter |7]], it is interesting to ask whether the NQA 
to the NJL model can be extended to finite T and /i. In this letter we shall demonstrate 
that the answer to this question is affirmative. We stress that the method presented here to 
extend the NQA to finite T and /i is quite general and is not restricted to the study of the 
NJL model. 

Our strategy is based on the extension of the thermal field theory formalism of Umezawa 
and collaborators || by including a finite chemical potential, and has been motivated by the 
application of the NQA to the BCS theory of superconductivity [[|. The basic idea is to apply 
a generalized thermal Bogoliubov transformation to the creation and annihilation operators 
of the asymptotic fields in the Haag expansion accompanied by a thermal doubling of the 
Hilbert space. In this way temperature and chemical potential are introduced simultaneously 
through the coefficients of the transformation and are treated as parameters on an equal 
footing. Quantities of interest are then obtained through the algebra of thermal creation 
and annihilation operators without specifiying the structure of the thermal vacuum state. 
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When the method is applied to the NJL model we recover the known mean field results 
for number and scalar densities, total chemical potential, gap equation and vacuum energy 
density using only the first term in the Haag expansion. Thus NQA in the lowest order Haag 
expansion is equivalent to the mean field approximation usually employed in the studies of 
the NJL model at finite T and \i. Furthermore, in the same approximation, we find that the 
Hamiltonian can not be diagonalized at finite T and in the broken symmetry phase of the 
model due to particle-hole and anti-particle and anti-hole excitations. This seems to be a 
general result valid for all versions of the NJL model at finite T and \i. 

In the following section we start by summarizing the pertinent results for the U(l) NJL 
model for (T, \x) = (0, 0) in the NQA. In Section III the NQA is extended to finite T and 
fi which is used to recover the known mean field results in Section IV. In the same section 
we present the off-diagonal Hamiltonian and show that it can only vanish in the T = 
or in the Breit-Wigner phase of the model. We repeat the same excercise in Section V 
using the SU(2) f ® SU(3) C NJL model including the scalar and axial vector color 3 diquark 
channels with their respective coupling constants. By demanding that the total Hamiltonian 
is diagonalized at T = a simple relation between these coupling constants is found. This 
relation is then used to show how the inclusion of the diquark channels can lower the vacuum 
energy density of the many-body system. In the final section we summarize the results and 
outline the procedure for going beyond the mean field approximation. 



II. NQA TO THE U(l) NJL MODEL: (7» = (0,0) 
We first consider the U(l) invariant version of the NJL model defined by the Lagrangian 



c = m>$m + 9o (^) 2 + (m^f 



(2-1) 



The charge conjugation symmetric form of the corresponding field equation of motion and 
Hamiltonian in momentum space are 



0>(q) = -y J dfadfadfrs b A (q-p x - Pi - Ps) 



[*(Pi),*(P2)]-,*(p 3 ) 



[*(Pl),75*(P2)]-,75*(P3 



(2.2) 



and 

H = (2tt) 3 

-(2tt) 3 J d%d%d%d% 5 3 (p 1 +p 2 + p 3 + p 4 )e-^ + ^ 



1 d% <S 3 (pi +p 2 )e-^ + ^^[*(pi),7 • (Pi -£)*(p 2 ) 



*(P1),*(P2) ^(P 3 ),*(P4) " *(pi),75*(P2) ^(P3),75^(P4) , (2.3) 



where ^f(p) = l I / ^(— p)7o, [A, B]± = AB ± BA and we shall adapt the notation d"p = 
d n p/(27r) 3 . When (T, /x) = (0,0) there are no contributions from the exchange (Fock) terms 
in the mean field approximation in this model. 
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For the model defined by Eq. ( |2.1| ) the Haag expansion for the Heisenberg field ^(p) is 
given by f|, 

+ E / S\q-Pi-P2)f {k) ( P uP2) : *m(Pi)Mn)^mfo)<W(Pa) : 

k J 

+ J d\>x<$p 2 5 4 (q~Pi - P2)e(pi,p 2 )C~ l : ^m(pi)S m (pi)E m (p 2 )S mk (p 2 ) : 

+ • • • . (2.4) 

Here \1/in is the asymptotic in-field of \l/ constrained to be on the mass-shell by the delta 
function 5 m {p) = <5 4 (p 2 — m 2 ). C = ij 2 lo is the charge conjugation operator and the 
symbol : : denotes normal ordering. The k summation runs over the ^275 \&, and 
^7^ bound states with corresponding on-shell in-fields Fjj$ and amplitudes f^(pi,p 2 )- 
Similarly e(pi,p 2 ) is the amplitude for the bound state with an on-shell asymptotic 
field Em. The ellipses represent higher order contributions involving terms with three or 
more products of ^m, ^in and asymptotic bound state fields together with corresponding 
amplitudes. An example of the third order contribution is given in ||. 

As shown in ||, inserting the first term in the above expansion into the operator equation 
of motion Eq. (|2.2|) yields the well known gap equation 



m 



4<7o / d 3 p — (2.5) 

J Up 



with L) p = \J\p\ 2 + m?. Note that the asymptotic mass m plays the role of the chiral gap. 
There are three solutions to this gap equation; m = corresponding to the phase with 
unbroken chiral symmetry and m = ±m for the broken symmetry phase. The vacuum 
energy densities for the two phases are 

_L (W) m=0 = _ 2% J d ^ d 3 q _ 2 JSp \p\ (2.6) 
1 {n) m=±m = _ 2gQ f d ^ d 3 q _ f d3p 2\p\ 2 + m 2 

V J J ujp 

To facilitate comparison with known results for finite T and \x we shall use the non- 
covariant regularization scheme and cut off the three momentum integral at \p \ = A. Then, 
as is well known, dynamical symmetry breaking occcurs only for those values of the coupling 
constant satisfying the relation 7c 2 /A 2 < go. For these values of go the vacuum energy density 
of the broken symmetry phase is lower than that of the unbroken phase, indicating that the 
phase with chiral condensates is the energetically favored phase. 



III. EXTENSION TO FINITE T AND fi 

In this work we shall explore the consequences of using only the first term in Eq. ( |2.4|) 
to describe the NJL model at finite T and \i. At first sight this seems to be a trivial 
approximation, but as we shall see it turns out not to be the case. Extension of the present 
work to include the second order terms will be discussed below. 
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The lowest order term in the Haag expansion for \1/ is [] 



m 



(2tt) 3 ^ 



Q s=± 



J2 b & s ) u (& s ) 6 (l° ~ u i) + rft (~^ s )v(-q, s)5(q° + u q 



(3.1) 



Because this term alone reproduces the gap equation Eq. (|2.5| ) and diagonalizes the Hamilto- 
nian, the operators b and d correspond to quasi-particle and quasi-anti-particle annihilation 
operators, respectively, which annihilate the interacting vacuum state at zero temperature 
and chemical potential |C?(0, 0)), 



b(q,s)\g(0,0)) = d(q,s)\g(0,0)) = 0. 



(3.2) 



Thus using the first order term in the Haag expansion is equivalent to subjecting the creation 
and annihilation operators in the Heisenberg field f to a Bogoliubov transformation into 
quasi-particle and quasi-anti-particle basis. 

However, the b and d operators do not annihilate the interacting vacuum state at finite 
T and \x denoted as \G(T, //)). In order to construct operators that annihilate \Q(T,fi)) we 
apply a generalized thermal Bogoliubov transformation to the b and d operators as well as to 
the annihilation operators corresponding to quasi-holes b and quasi-anti-holes d as follows 



b(q, s) = a q B(q, s) - s(3 q B^(-q, s) 
b(q, s) = a q B(q, s) + s/3 q B\-q, s) 
d(q, s) = >y q D(q, s) - s5 q D ] (-q, s) 
d(q, s) = i q D(q } s) + s5 q D\-q, s) 



(3.3) 
(3.4) 
(3.5) 
(3.6) 



Here B and B annihilate a quasi-particle and a quasi-hole at finite T and /i, respectively, 
and D and D are the thermal annihilation operators for quasi-anti-particles and quasi- 
anti-holes, respectively. These operators annihilate the interacting thermal vacuum state 
for each T and fi. 

B(q,s)\g(T, f i)) = B(q,s)\g(T, f x))=D(q,s)\g(T^)) = D^ (3.7) 

In the NQA it is not necessary to specify the structure of \Q(T,[i)). We only have to 
assume that it is annihilated by the annihilation operators in the Haag expansion. The 
thermal doubling of the Hilbert space accompanying the thermal Bogoliubov transformation 
is implicit in Eq. ( |3.7D where a ground state which is annihlated by thermal operators B, 
B, D and D is defined. 

In addition, thermal operators satisfy the Fermion anti-commutation relations 



(27r)^6 3 (p-q)6 Slt 
m 



B\p,s{),B{q,Si) 
B\p, Sl ),B{q,s 2 ) 



D\p, Sl ),D(q,s 2 ) 
D\p, Sl ),D(q,s 2 ) 



(3.8) 



1 We use the conventions and normalizations of Itzykson and Zuber [1C|. 
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with vanishing anti-commutators for the remaining combinations. The coefficients of the 
transformation are a q = Jl — n~, (3 q = <Jnq~, 7q = \Jl — nf and 5 q = yjnf, where = 



,K±M)/(fc s T) 



+ 1] 1 are the Fermi distribution functions for particles and anti-particles. 



They are chosen so that the total particle number densities are given by 



1 m 



-^(g(T,ti)\b\q,s)b(q,s)\g(T,ri) 



1 m 



--(g(T,ti)\<#(q,s)d(q,s)\G(T,ri) 



Vu, 



(3.9) 
(3.10) 



Hence in this approach temperature and chemical potential are introduced through the 
coefficients of the thermal Bogoliubov transformation and are treated as parameters. 
The first order Haag expansion for \l/ extended to finite T and fi is 



m 



q s=± 



Y,\ a q B(q,s) - s/3 p B\-q,s) u(q, s)5(q° - u q ) 



+ 



lq D\-q,s)-sS q D(q,s) v(-q, s)5(q° + u q ) . (3.11) 



Note that this extension is independent of the structure of the Lagrangian and is therefore 
not restricted to the NJL model. If bosonic fields are present, as in the second order terms 
in Eq. fl2.4|) , a thermal Bogoliubov transformation for bosonic creation and annihilation 
operators must be introduced together with a redefinition of the thermal ground state so 
that it is annihilated by both fermionic and bosonic annihilation operators. In the following 
sections we shall use the above ansatz for the ^ field to study the NJL model at finite T 
and [i. 



IV. NQA TO THE U(l) NJL MODEL: (T,fi) + (0,0) 

Before extending the U(l) NJL model to finite T and /i we first calculate the number 
(AT) and scalar [S] densities using the ansatz given in Eq. (|3.11|) . It is a straightforward 
exercise to show that the number density is given by 



V 

2 / d 3 p 



-5 2 

V V 



(4.1) 



while the scalar density is found to be 
1 



S = -{2-nY / d%d% &\p 1+ mG{T^)\- *( Pl ),*(p 2 ) \G(T,fi)) 



-2jd?p^{l-Pl-5f) 



(4.2) 



We now insert a bare chemical potential term yUo^7o^ in the Lagrangian Eq. ( [2.1|) and 
use Eq. ( |3.11|) in the corresponding equation of motion and renormal order. Keeping only 
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the linear terms in the thermal field operators, one finds that the equation of motion is given 
by 



(m + /i o7 o)^(g) = I [4g, 



dp — (1 — (3, 



(4.3) 



Whereas the exchange terms did not contribute at (T, //) = (0,0), when introducing T and 
\i through the thermal Bogoliubov transformation one finds contributions to the equation 
of motion from both the direct and exchange terms. In fact, the exchange terms generate a 
contribution to the chemical potential proportional to the number density. From Eq. ( |4.3| ) 
we see that the total chemical potential is given by /i = /i — g§N ' . 

The gap equation for (T, /i) ^ (0,0) is found by equating the scalar terms in Eq. 
The results is 



m 



4so/d 3 P^(l-/3 p 2 -5 ; 



(4.4) 



and reduces to Eq. ( |2.5| ) in the (T, /i) — > (0,0) limit. To obtain the vacuum energy density 
e(T,fi) we use Eq. (|2.3| ) (plus the /i term) and keep terms with two contractions after 
renormal ordering. We find 



e(T,ii) = ±{g(T,ii)\H\g(T,ri) 



2g J d 3 pd 3 q + ^g M 2 + 2 J d 3 p cu p (1 - (3 2 p - 5 2 p ) - ^- + \ 



m 



(4.5) 



Above expressions for the number and scalar densities, total chemical potential, the gap 
equation and the vacuum energy density all agree with the standard mean field results for 
the U(l) NJL model 0]. We note that unlike in the variational method the explicit structure 
of the thermal vacuum state need not be specified when obtaining these results. Instead 
they were determined by the algebra of the thermal fields defined by Eq. ( |3.7|) and ( |3.8D . 
We now turn to the discussion of the diagonalization of the model Hamiltonian. 

In order to obtain the off-diagonal terms in the Hamiltonian, denoted as Hod, to sec- 
ond order in the thermal creation operators we keep terms involving one contraction when 
renormal ordering Eq. ( |2.3| ) (again, with the /xq term), and exploit the fact that the ther- 
mal annihilation operators annihilate the thermal vacuum state. Using the gap equation 
Eq. ( f4.4| ) the result can be written as 



Hod = m^2sjd 3 p^ a p (3 p B\p, s)B\-p, s) + lp 5 p D\-p, s)D\p, s) 

a p (3 p B\p, s)B\-p, s) - lp 8 p D\-p, s)D^p, s) 



(4.6) 



We see that the off-diagonal terms consist of particle-hole and anti-particle-anti-hole exci- 
tations coupled to the spin triplet state. These terms vanish in the chirally symmetric phase 
when m = 0, or in the zero temperature limit when a p /3 p = 7 P 5 P = 0. However Eq. ( ^4.6| ) can 
not vanish for finite T in the broken symmetry phase. 
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In the standard treatment of BCS theory of superconductivity |TT[ the gap equation 



is obtained by demanding that the second order off-diagonal Hamiltonian vanishes. We 
could also have followed this path and derived the gap equation from Hob- Instead we used 
the result obtained from the equation of motion to simplify the second order off-diagonal 
Hamiltonian as a consistency check. Also, in the BCS theory the condensation takes place 
between particles and holes in the spin singlet state and thus contributions from spin triplet 
terms in the Hamiltonian are usually ignored JTJJ . In the present case the chiral condensate 
is also a spin singlet bound state and if we adapt an analogous approximation Eq. ( |4.6| ) can 
also be ignored. 



V. SU(2) F ® SU(3) C NJL MODEL AT FINITE T AND y, 

We now apply the method presented in the preceeding sections to the SU(2) f <g> SU(3) C 
NJL model consisting of scalar (0 + ,T = 0) and pseudoscalar (0 + ,T = 1) color singlet qq 
channels as well as scalar and axial vector (1 + ,T = 1) color 3 qq channels. This model was 
used in [O] to study three quark bound states in the NJL model and is defined as 



+ 92- 



(*7 5 Cr 2 X A V T ) (* T C- 1 7 5-r 2 A A *) 



+ 93 2 



3r = 



(5.1) 



where X A with A = 2,5,7 project out the color 3 channel.0 Summation convention is 
implied for repeated indices. Positive values for the coupling constants gx, g 2 and g^ implies 
attraction in color 1 and 3 channels. 

Just as in the case of the U(l) NJL model we Haag expand the Heisenberg field operators 
and apply the thermal Bogoliubov transformation to the first order terms. Use of the 
equation of motion corresponding to Eq. (|5.1|) immediately yields the total chemical potential 
and the gap equation for finite T and \i 



71 = 26*?! / d^^{l-(3 2 -5 2 p 



(5.2) 
(5.3) 



These results for fi and the gap equation are identical to the Hartree-Fock results for the 
SU(2)f (g) SU(3) C NJL model without the diquark channels |7|]. Indeed, we find that to the 
lowest order in the Haag expansion the diquark channels do not contribute to the equation 
of motion. 

The diquark channels do contribute to the vacuum energy density which, with the use 
of Eq. (|5.2| ) and ( |5.3|) , can be written as 



2 The Gell-Mann matrices A a are normalized as Tr(A a A b ) = 26 a t,. 
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e(T, = -6{ I2 9l I d 3 pd 3 q + 2 / d 3 p co p (1 - (3 2 p - 5 2 p ) - (26,! + 6, 2 - 72g 3) 



ni 



26 9l 



+ - (6, 2 - 36, 3 ) 

4,i 



/ d 3 p - 51 



(5.4) 



Above expression reduces to the known vacuum energy density for the SU(2)f (g) SU(3) C NJL 
model when g 2 — ,3 — 0. Similarly, the second order off-diagonal Hamiltonian can be 
simplified to the following form 



T^od = m / d 3 p <J2 S 



I 777 

1 + ^— (6, 2 + 72,3)^ 



26, 



a P P P B L(P, s)Bl(-p, s) + j p S p Dj e (-p, s)Dj e (p, s) 



-5>- 



s Up i 



H + (6, 2 + 36, 3 ) J d 3 p 



^ ( %-?p 



E 

SI ,52 



j fjl^ 
(6,2 + 72,3) — 



26, 



UJ, 



p j 



a P 7pD}M s 2 )Bl(-p, si) - s^PpSpDl^-p, s x )B} e (p, s 2 ) 



(5.5) 



where I = u,d and e = 1, 2, 3 are the flavor and color indices, respecitvely and Xi is the two 
component spinor for spin Sj. 

One sees immediately that the Hamiltonian for the SU(2)f ® SU(3) C NJL model is exactly 
diagonalized in the Breit-Wigner phase when m = 0. However in the broken symmetry 
phase the first two terms in Eq. ( p.5[ ), which is similar to Eq. Q4.6| ), vanish only in the zero 
temperature limit but the third term in general does not. This term consists of particle- 
anti-particle and hole-anti-hole excitations with zero total three momentum with the spins 
coupled to both and 1. Since it contains spin singlet excitations it is desirable that this 
term vanishes exactly in order to obtain a true ground state of the many-body system at 
least in the zero temperature limit. This can be accomplished if we impose the condition 
,3 = —3^,2 on the coupling constants for the scalar and axial vector diquark interactions. 
An immediate consequence of this restriction is that the scalar and axial vector qq channels 
can not both be attractive or repulsive. 

In Figure la we plot the phase transition curve for the SU(2)f ® SU(3) C NJL model 
using the Hartree approximation for the qq channel .0 The input parameters for all our 



3 In this approximation Eq. ( |5.2| ) becomes \x = \iq and the coefficient infront of the integral in 
Eq. (|5.3|) changes from 26 to 24. 
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numerical work are A = 0.05 GeV, g\ = 5.01 GeV 2 and g?. = —I2g^ = 3.11 GeV 



which are taken from |T3[]. On the phase transition line m = and the model Hamiltonian 
is exactly diagonalized. The corresponding vacuum energy densities, e(T, /i) — e(0,0), are 
shown in Figure lb. The upper curve corresponds to the energy density obtained with only 
the qq channel while the lower one includes both the qq and qq channels. For all values of 
temperature and chemical potential on the phase transition line the inclusion of the scalar 
and axial vector diquark channel lowers the vacuum energy density of the system. 

The effect of the diquark channel on the energy density is similar in the broken symmetry 
phase at finite temperature. Figure 2a shows the behaviour of the asymptotic mass m as a 
function of fi for T = 0.1 GeV. The phase transition is second order, which is a consequence 
of the Hartree approximation for the qq channel, and the critical chemical potential is found 
to be yU c = 0.253 GeV. Below fi c the Hamiltonian is only approximately diagonalized and in 
Figure 2b we plot the vacuum energy density as a function of fi. As in Figure lb the upper 
curve is obtained with only the qq channel while the lower curve has contributions from both 
the qq and qq channels. Note that in both cases the values of the vacuum energy densities 
are largest at fi = fi c indicating that the broken symmetry phase is the energetically favored 
phase. We also find that the addition of the qq channel lowers the vacuum energy density 
for finite density at zero temperature. Hence the inclusion of the diquark channel can lower 
the vacuum energy density of the system when the Hartree approximation is invoked for the 
qq channel. 



VI. SUMMARY AND OUTLOOK 

In this work we have extended the NQA to quantum field theory to finite temperature 
and chemical potential. The basic idea is to subject the asymptotic fields in the Haag ex- 
pansion to a thermal Bogoliubov transformation accompanied by a thermal doubling of the 
Hilbert space. Temperature and chemical potential are introduced non-linearly through the 
coefficients of the transformation. The NQA allows one to solve the operator equations of 
motion without explicitly specifying the structure of the ground state although it can be 
obtained a posteriori by diagonalizing the Hamiltonian. All quantities of interest are deter- 
mined through the algebra of asymptotic fields. One only needs to assume the existence of a 
ground state which is annihilated by the annihilation operators used in the Haag expansion. 
The thermal doubling of the Hilbert space is implicit in Eq. ( |3.7|) where a thermal ground 
state was introduced which is annihilated by the thermal annihilation operators. 

We tested our formalism on the U(l) and SU(2) f <g> SU(3) C versions of the NJL model 
and found that the first term in the Hagg expansion can reproduce the mean field results for 
the number and scalar densities, the chemical potential, the gap equation and the vacuum 
energy density for finite T and fi. Furthermore we found that the second order Hamiltonian 
can not be diagonalized at finite T when the models are in the broken symmetry phase. 

When scalar and axial vector qq channels are added to the SU(2) f ® SU(3) C NJL model 
we found it necessary to impose a relation on the scalar and axial vector coupling constants 
to diagonlize the Hamiltonian for zero temperature. As a consequence scalar and vector 
diquark channels can not simulateously be attractive or repulsive. For an attractive scalar 
diquark channel and in the Hartree approximation to the qq channel, this restriction can lead 
to the lowering of the vacuum energy density. These findings should be taken into account 
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when the work developed in |12] to study three quark bound states in the NJL model is 
extended to finite T and //. 

To go beyond the mean field approximation it is necessary to add higher order terms 
in the Haag expansion as shown in Eq. ( p.4|) . Here the second order terms involve both 
fermionic and bosonic fields, the latter corresponding to various bound states in the model. 
It is necessary to subject both types of fields to their respective thermal Bogoliubov transfor- 
mations accompanied by a redefinition of the thermal ground state so that it is annihilated 
by both the fermionic and bosonic annihilation operators. In the NJL model it is possible 
to analytically calculate the bound state amplitudes by exploiting their behaviour under 
parity, time-reversal and charge-conjugation transformations ||. Once these amplitudes 
have been determined for finite T and /i they can be fed back into the Haag expansion of 
the Heisenberg fields and quantities of interest recalculated. 

We emphasize that the formalism presented here is quite general and has wide applica- 
tions. Recently the importance of the scalar diquark channel in the study of color supercon- 
ductivity has been emphasized by several authors using an effective four-Fermion instanton 
induced interaction JR]. Our method is very well suited for this application and investiga- 
tion into the phenomenology of color superconductivity at finite temperature and chemical 
potential beyond the mean field approximation is currently in progress. Other applications 
include studies of pion or kaon condensations or properties of pion gas in chiral perturbation 
theory. In fact, our method can be applied to any quantum field theory in Minkowski space 
with well defined asymptotic fields. 
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FIG. 1. a) Phase trasnsition line for the SU(2) f <g> SU(3) C NJL model with A = 0.05 GeV and 
<7i = 5.01 GeV -2 . b) Vacuum energy density e(T, /z) — e(0, 0) on the phase transition line begining 
at (T,n) = (0,0.307) GeV and ending at (T,fi) = (0.185,0) GeV. The upper and lower curves 
have been obtained with g\ = 5.01 GeV -2 , Q2 = <?3 = and g\ = 5.01 GeV -2 , §2 = — I2gs = 3.11 
GeV -2 , respectively. 
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FIG. 2. a) Asymptotic mass m in the SU(2)f (g> SU(3) C NJL model as a function of chemical 
potential fi for T = 0.1 GeV with the same model parameters as in Figure la. b) Vacuum energy 
density e(T, fi) — e(0, 0) as a function of /i for T = 0.1 GeV. The upper and lower curves have been 
obtained with gi = 5.01 GeV~ 2 , g 2 = g 3 = and g x = 5.01 GeV" 2 , g 2 = -12g 3 = 3.11 GeV" 2 , 
respectively. 
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